Phase Transitions in Nucleonic Matter and Neutron-Star Cooling 
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A new scenario for neutron-star cooling is proposed, based on the correspondence between pion 
condensation, occurring in neutron matter due to critical spin-isospin fluctuations, and the metal- 
insulator phase transition in a two-dimensional electron gas. Beyond the threshold density for pion 
condensation, where neutron-star matter loses its spatial homogeneity, the neutron single-particle 
spectrum acquires an insulating gap that quenches neutron contributions to neutrino-production 
reactions and to the star's specific heat. In the liquid phase at densities below the transition point, 
spin-isospin fiuctuations are found to play dual roles. On the one hand, they lead to a multi-sheeted 
neutron Fermi surface that extends to low momenta, thereby activating the normally forbidden 
direct-Urca cooling mechanism; on the other, they amplify the nodeless P-wave neutron superfiuid 
gap while suppressing S-wave pairing. In this picture, lighter stars without a pion-condensed core 
experience slow cooling, while enhanced cooling occurs in heavier stars through direct-Urca emission 
from a narrow shell of the interior. 
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Observations designed to detect and measure thermal 
radiation from neutron stars with ages between 10^ and 
10^ years can provide valuable information on the inter- 
nal temperatures T of these stars at the stage when cool- 
ing proceeds via neutrino emission from their neutrino- 
transparent cores 1]. Such information may be used to 
constrain theories of hadronic matter at high density. 
The relevant neutrino reactions are (i) the so-called direct 
Urea (DU) process, in which beta decay and inverse beta 
decay operate in tandem on thermally activated nucle- 
ons and electrons, (ii) the modified Urea (MU) process, 
involving a nucleon spectator in the Urea reactions, (iii) 
neutrino bremsstrahlung from nucleon- nucleon collisions, 
and (iv) neutrino emission due to Cooper pairing in the 
superfiuid phase Of these, only the DU mechanism 

can produce rapid cooling, as the data seems to require 
for some stars. However, its operation is normally forbid- 
den by the large mismatch of neutron and proton Fermi 
momenta, which can be compensated by spectators in 
the MU process generally associated with slow cooling. 
The importance of the Cooper-pairing mechanism (iv) is 
highly sensitive to the proton and neutron gaps, whose 
values in turn depend crucially on the internal composi- 
tion and structure of the star. 

A new perspective on neutron-star interiors is sug- 
gested by laboratory studies of (i) the metal-insulator 
transition in two-dimensional (2D) silicon samples of low 
disorder J4, and (ii) the liquid-solid phase transition in 
2D ^He |5|- Quite far from the transition point, the elec- 
tron and ^He systems involved may be described as nor- 
mal Fermi liquids obeying Landau theory. Beyond the 
transition, these systems become inhomogeneous, with 
an insulating gap in the single-particle (sp) spectrum. 



Similar behavior of the sp spectra has been observed in 
the "normal" states of high- Tc superconductors 6]. Cold 
neutron matter also behaves as a normal (or superfiuid) 
Fermi liquid at relatively low densities. With increas- 
ing depth in the star and coincident increase in baryon 
density p, spin-isospin fluctuations grow in importance. 
At a critical density pc, a condensate of spin-isospin ex- 
citations, specified by critical wave number qc ^_PEjii 
forms in the channel with tt^ quantum numbers |3, 
State-of-the-art microscopic calculations HI predict that 
tt" condensation (PC) sets in at a density comparable to 
the equilibrium density po = 0.16 fm~^ of symmetrical 
nuclear matter, while earlier phenomenological estimates 
put pc in the range 0.3 - 0.4 fm~^ 0jl3 • Comparing these 
values to the central density p(0) ~ (0.5 — 1.0) fm~^ 
of a typical neutron star, we infer that at low T, a signif- 
icant part of the stellar bulk exists in the inhomogeneous 
PC phase, rather than as a Fermi liquid. 

Accordingly, in contrast to the pionic cooling scenario 
reviewed in Ref. [i(tI | , we propose that in the domain oc- 
cupied by the PC, the neutron spectrum acquires an in- 
sulating gap exceeding in value any expected superfiuid 
gaps in neutron matter. It follows that in this region of 
the star all cooling processes requiring the participation 
of neutrons are strongly inhibited. The neutron contribu- 
tion to the specific heat C{T) is suppressed as well, since 
its leading term ~ now derives from the phonon spec- 
trum rather than from the gapped sp excitations. The 
properties of the embedded proton subsystem are not 
affected so dramatically. Analogously to the electronic 
subsystem in solids, the protons form a "conductivity 
band" with a sp spectrum of unaltered shape, specified 
by an effective mass, and the behavior of their contribu- 
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tion Cp(T) ~ T to the specific heat remains unchanged. 
Thus, in the new coohng scenario, the large region of 
the star in which the tt" condensate holds sway becomes 
irrelevant to the cooling process, except for the effects 
of neutrino-generating reactions involving protons as the 
only nucleonic participants, along with the proton con- 
tribution to the stellar specific heat. 

As the internal temperature drops, the proton subsys- 
tem undergoes a superconducting phase transition, which 
we assume to originate from phonon-induced attraction. 
A rough estimate of the critical temperature, ~ 20 
keV, is then obtained in terms of the proton Fermi en- 
ergy in the same way as for ordinary superconductors. 

Having addressed the organization of neutron-star 
matter at moderately high densities, let us now analyze 
what is transpiring in the domain 0.5 po < p < pc cor- 
responding to the outer core, where the stellar material 
is a neutron liquid with fluid admixtures of protons and 
neutralizing leptons. We argue that in this region, spin- 
isospin fluctuations have a strong impact on neutron pair- 
ing, suppressing it in the S'-wave channel and enhancing 
it in P-waves. This situation is familiar from the physics 
of superfluid ^He, where spin fluctuations play the key 
role in promoting P-pairing over ^-pairing [llj. To esti- 
mate the neutron triplet-P gap, we adopt the BCS for- 
malism and write the gap equation as 



A(p) = - y [V(p, p')+V- (P-P'; oj = E{p'))] ^^^dr' . 

(1) 

with dr = d?p/ {2tt)'^. In the triplet P-wave channel, the 
gap function has the form A(p) = idikPkO'ia2, with coeffi- 
cients dik yet to be determined. The quasiparticle energy 

E{p) is given by E{p) = [^^(p) -f Tr(A(p)At(p))] 

where ^{p) is the sp spectrum of the normal Fermi 
liquid relative to the chemical potential /z, and 
Tr(A(p)A+(p)) = + E™«m>^2m(n). In triplet P- 
pairing, the regular interaction V is nearly exhausted by 
the frequency- and density-independent spin-orbit com- 
ponent of the scattering amplitude. The fluctuation 
contribution to the pairing interaction is written as 

Kf3,-ys(^^^) = ^«(9)('^a7 ' q) Re D{q, u;){<T ps ■ q)/ml, 
where A„ is an effective charge accounting for the renor- 
malization of the associated vertex part. The propa- 
gator D is given by the Migdal formula —D^^{q,uj) = 
+ ml + Uisri{q,u} = 0) + Ui^]y{q,u!), made up of the 
ordinary part njvAr(q,w) = —Bq^ — iq\uj\M'^ / {2'!Tm1) of 
the pion polarization operator 11, along with the term 
njv/(g,a; = 0) = —q^pj pi[\ -f q^ lq\) arising from pion 
conversion into a A-isobar and neutron hole. In the do- 
main of critical fluctuations, we have 01 
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where rj = {pc — p)/pc- The parameters 7 and k are de- 
termined from the obvious relations (1 — B)q'^ + — 



Qc'^cCc = and 1 - B - r^Q = 0- with = Pc/pi 
and (c — (1 + 'Zc/'?/)^^- Simple algebra leads to 
qc - (m,(7/)'/' (1 - B)-^/\ = {m^/qi + VT^)\ 
7^ — TcCci and — VcCcqc/Qi- Employing the parame- 
ter set q] = 5ml., pj ~ 1.8/90, and B ~ 0.7 from Ref. 0, 
one arrives at qc ~ 1.9™^, pc — 2 po, 7 ~ 0.4, and 
K ~ 0.7. As will be seen, the contribution to the gap 
value from spin-orbit forces is insignificant; with its ne- 
glect Eq. ^ is finally recast as 



dtkPk = - J 



A?. (9) 
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where q = p — p'. In arriving at this result, we 
have made use of the relations 'Ylif,k'^i'^i'^'i)aii{,'^rn)!i& — 

Unfortunately, the pairing problem ^ still defies full 
solution. In neutron matter, only the spectrum of so- 
lutions restricted to the ^P2-^P2 channel has been ex- 
plicated 01 . The solution set contains both nodeless 
and nodal combinations of the different basis states. At 
T <C Pc, a nodeless solution wins the energetic compe- 
tition, and we anticipate that this feature, also inherent 
in superfluid ^He, is present here as well. For nodeless 
solutions, an adequate approximation to the gap in the 
sp spectrum can be obtained by retaining only the A^ 
term in Tr(A(p)A"'"(p)). The matrix dik then becomes 
proportional to 5ik and the angle integration is obviated, 
giving rise to the expression q^ (p • p') — 2(q • p)(q • p') 
in the numerator of the r.h.s. of Eq. (|3Jl. Furthermore, 
the system (3) becomes decoupled, and we are left with a 
single integral equation to solve. Exploiting the fact that 
the propagator D{q,uj = 0) is peaked at q = qc < 2pi?„, 
this equation is simplified to 
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where = X^qc). In deriving this formula, we 
have made the replacements p • p' = {p"^ + (p')^ ^ 
q2)/2->p|,-92/2, q • p = p2 _ p . p'^g2/2, and q • p' = 
p • p' — (p')^^ — ql/2; their validity has been confirmed in 
numerical calculations. Since the parameter A„ is as yet 
uncertain, numerical calculations of the gap A(p) have 
been performed for three different values, with the re- 
sults shown in Fig. ^ We see that the P-wave neutron 
gap is dramatically magnified in comparison with stan- 
dard estimates of 5-wave and P-wave gaps that ig- 
nore spin-isospin fluctuations, justifying our neglect of 
the regular interaction. 

Next we assess the inhibitory effect of spin-isospin 
fluctuations on proton S'-wave pairing. To estimate 
the suppression factor, we appeal to the property 
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FIG. 1: Neutron P-wave gap A„ at the Fermi surface (in 
MeV) (panel (a)) and suppression factor Ap/Ap"' for proton 
S-wave gap (panel (b)), versus baryon density p in units of 
pc- (The proton fraction is taken as 0.06.) Values of and 
Ap are indicated by numbers near the corresponding curves. 



{S=Q\<j](j1\S^Q) = ~5ik and multiply both sides of 
Eq. J^l, rewritten for the full proton S'-wave gap function 

1/2 



Ap, by A^,"V2i^ob), where Eo{p)= [e{p) 

and Ap°^ is the gap value found 0| when critical fluc- 
tuations are neglected. Integration over the intermedi- 
ate proton momentum and simple manipulations utiliz- 



ing the gap equation for Ap 
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with p = p(^) determined from the formula for the sp 
spectrum ^(p). Results of numerical calculations of the 
ratio A/A*^°) based on Eq. ||SJ) are also presented in Fig.D 
Proton S-pairing in the liquid phase is suppressed over a 
wide density range. 

The above results suggest that neutron P-pairing in 
the liquid domain of the core is enhanced so strongly 
by spin-isospin fluctuations that the neutrino emissiv- 
ity due to neutron Cooper pairing, which behaves as Q 
QnCp ~ i^/T)^ exp{—2A/T), is completely suppressed 
by the exponential factor. At the same time, recent mi- 
croscopic calculations indicate that the ^Sq proton 



gap in neutron-star matter depends crucially on the den- 
sity p, falling off rapidly when p exceeds po- In view 
of the suppression factor from spin-isospin fluctuations 
found here and plotted in Fig. 1, we conclude that pro- 
ton pairing is irrelevant to the neutrino-cooling stage. 

We now turn to the role of the direct Urea process. 
Despite their limitations, the available experimental data 
on the surface temperatures Ts of neutron stars give evi- 
dence for the existence of slow and rapid cooling tracks. 
It is generally presumed that DU reactions are somehow 
involved in the rapid cooling process. Yet if one adopts 
the best available equations of state of neutron matter 
derived from first principles, this highly efficient cool- 
ing mechanism is precluded in all but the most massive 
neutron stars, because of the large difference between 
neutron and proton Fermi momenta. However, a novel 
route to the DU process in dense matter is opened by a 
rearrangement of the neutron quasiparticle distribution 
n(p) at a critical density pr < Pc- The rearrangement is 
precipitated by critical spin-isospin fluctuations seething 
in the neutron liquid near its inner boundary with the 
PC domain Ts'l. 

To elucidate this phenomenon, we employ the Landau 
relation for the sp spectrum e{p) in the specific form 



de{p) deoip) 



dp 
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V-(p-pr,. = 0)^dn, 
dpi 
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where eo{p) — p^ /2M^ is the regular part of the neutron 
spectrum with A/q ~ O.TAf, the customary value of the 
neutron effective mass in the absence of spin-isospin fluc- 
tuations. Upon straightforward momentum integration, 
relation yields a closed RPA-like equation 
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well suited to investigation of the rearrangement of the 
Landau state as the density climbs to the critical value pc- 
Results of numerical calculations based on Eq. (jS]), de- 
picted in Fig. 2, indicate that the neutron Fermi surface 
becomes doubly connected at rjr — {pc — Pr)/ Pc — 0.065. 

With the redistribution of quasiparticles in momentum 
space, an inner neutron Fermi surface emerges at a low 
momentum pi, enabling momentum/energy conservation 
in the DU reactions and unleashing rapid cooling. When 
the neutron liquid becomes supcrfluid under decreasing 
temperature, a triplet pairing gap A.^ forms on the in- 
ner Fermi surface and applies the brakes to DU emission. 
Even so, we find that this novel opportunity for DU cool- 
ing 12 is not destroyed. To gauge the neutron gap value 
Ai at the inner Fermi surface, we again employ Eq. 
Numerical computation yields 



Mpf)Pi/pf 



(9) 
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FIG. 2: Neutron quasiparticle spectrum ^(p) in units of 
el = pI/2M = (37rVc)^''V2M, plotted for different densi- 
ties. Numbers near the curves give the corresponding values 
of = (/5c — p)/pc. Inset shows the quasiparticle spectrum 
^(p) evaluated for 77 = 0.062, together with the associated 
quasiparticle momentum distribution n(j>). 

The DU process may operate with full force only if the 
gap value A„ is markedly less than 100 keV. According to 
Eq. ((SJ, this constraint is met at the inner Fermi surface 
if Pi{p) < 0.02pi?. The numerical calculation underlying 
Fig. 121 indicates that this inequality does hold in a thin 
shell of stellar material where p — pr ^ 5 x 10~^/9c. 

Now consider, in outline, how our model explains ex- 
perimental cooling data. We assume that neutron stars 
such as RX J0822-4300 and PSR B1055-52 have rel- 
atively low masses, such that pion condensation is not 
triggered in their cores. We then attribute the high val- 
ues of their surface temperatures to the absence of cool- 
ing mechanisms other than neutrino bremsstrahlung from 
pp collisions, emission processes involving neutrons being 
switched off by substantial neutron P-pairing. 

In explanation of the enhanced cooling rates of the 
Vela, Geminga, and 3C58 pulsars, we suggest that their 
central densities exceed some 0.5 fm^^, and hence that a 
pion condensate occupies a significant portion of their 
interiors. Rearrangement of the neutron quasiparticle 
distribution in a small region adjacent to the boundary 
between liquid and solid phases lifts the ban on the DU 
process. However, the DU cooling rate, proportional to 
the volume where this mechanism operates vigorously, 
is suppressed due to the restriction of the process to a 
narrow shell. Consequently, the corresponding values 
need not lie far below those of the first group of neutron 
stars. One feature of our model that warrants further 
examination is the fast transition from slow to enhanced 
cooling under increasing stellar mass. 

Restriction of DU reactions to a narrow shell does not 
apply for the most massive neutron stars, in which the 



internal pressure becomes high enough to melt the pion- 
condensate lattice. In the dense liquid core, the DU re- 
actions are allowed to proceed apace even if the standard 
microscopic equation of state |2| is employed. As yet, 
no observational evidence exists for such an unrestrained 
DU cooling mechanism. 

In this letter, we have proposed a new picture of the 
interior of neutron stars based on the similarity between 
pion condensation and the metal-insulator phase transi- 
tion in the 2D strongly-correlated Fermi systems. We 
have demonstrated that the incorporation of spin-isospin 
fluctuations has dramatic effects on the neutron and pro- 
ton gaps in the liquid part of the stellar interior and in- 
duces a rearrangement of the neutron Fermi surface that 
triggers the direct-Urca reaction. The associated cooling 
scenario provides for slow and accelerated cooling tracks 
that do not conflict with observational data. 
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